Ratchet universality in the presence of thermal noise 
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We show that directed ratchet transport of a driven overdamped Brownian particle subjected to 
a spatially periodic and symmetric potential can be reliably controlled by tailoring a biharmonic 
temporal force, in coherence with the degree-of-symmetry-breaking mechanism. We demonstrate 
that the effect of finite temperature on the purely deterministic ratchet scenario can be understood 
as an effective noise-induced change of the potential barrier which is in turn controlled by the 
degree-of-symmetry-breaking mechanism. A possible experimental realization of the present ratchet 
scenario is proposed in the context of a nanoscale ratchet composed of a pinion and a rack coupled 
via the lateral Casimir force. 
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Directed transport without any net external force, the 
ratchet effect [1,2], has been an intensely studied inter- 
disciplinary subject over the last few decades owing to 
its relevance in biology where ratchet mechanisms are 
found to underlie the working principles of molecular mo- 
tors [3,4], and to its wide range of potential technological 
applications including micro- and nano-technologies. Di- 
rected ratchet transport (DRT) is today understood to be 
a result of the interplay of symmetry breaking [5] , nonlin- 
earity, and non-equilibrium fluctuations, in which these 
fluctuations include temporal noise [2], spatial disorder 
[6], and quenched temporal disorder [7]. In extremely 
small systems, including many of those occurring in bio- 
logical and liquid environments as well as many nanoscale 
devices, DRT is often suitably described by overdamped 
ratchets, in which inertial effects are negligible in com- 
parison with friction effects. Three of the great diversity 
of contexts in which overdamped ratchets have been con- 
sidered are annular Josephson junctions embedded in an 
inhomogeneous magnetic field [8], synchronization phe- 
nomena of coupled oscillators [9], and quasiperiodicity 
routes in cold atoms [10]. Recently, the dynamics of 
nanoscale systems composed of one corrugated cylinder 
(pinion) and one [11] or two [12] corrugated plates (racks) 
at zero temperature have been studied in the overdamped 
regime. The pinion and rack(s) have no mechanical 
contact, but are coupled via the (lateral) Casimir force 
[13]. Indeed, it has been shown that the lateral Casimir 
force between corrugated surfaces provides a possibility 
for frictionless transduction of lateral forces in nanome- 
chanical devices without any physical contact between 
them [14,15]. More recently, there have been studies of 
a Casimir-force-driven ratchet in the presence of inertia 
and finite thermal noise [16], and of a thermal Brownian 
motor coupled by Casimir interaction in the overdamped 
regime [17]. 




FIG. 1: (Color online) The schematics of the noncontact rack 
and pinion, with the rack vibrating biharmonically laterally. 
The rack and pinion have sinusoidal corrugations of the same 
amplitude and wavelength A, while r and R are the internal 
and external radii of the pinion. 



by studying a strongly damped pinion at finite temper- 
ature that is kept at a distance from a vibrating rack 
(see Fig. 1). For deterministic ratchets, this has been 
shown to also be the case for topological solitons [7] and 
matter- wave solitons [20] . We consider the standard case 
in which both the pinion and the rack present sinusoidal 
corrugations with a single wavelength A, and assume, in 
contrast with the previously considered cases of uniform 
[21] and harmonic [22] motion, a biharmonic motion of 
the rack: 



y = A[ri cos (ujt) + (1 — 77) cos {2Ldt + ip)] 



(1) 



where A is an amplitude factor, and the parameters 77 g 
[0, 1] and if account for the relative amplitude and initial 
phase difference of the two harmonics, respectively. The 
rack and pinion are coupled by the lateral Casimir force 
between the neighbouring surface areas. 



In this Letter, we show how ratchet universality [18,19] 
works subtly in the context of noisy overdamped ratchets 
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where F is the amplitude and x — y is the lateral relative 
displacement. After using x = R9, where 9 is the angle 
of rotation, one obtains the equation of motion for the 
pinion: 



RF sin 



27r ■ 



= (3) 



where C is the rotational friction coefficient, ^ (t) 
is a Gaussian white noise with zero mean and 
{Ht)^{t + s)) = S{s), and D = 2ChT with h and T 
being the Boltzmann constant and temperature, respec- 
tively. After substituting Eq. (1) into Eq. (3) and using 
the dimensionless variables 
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one can rewrite the equation of motion as 



(r) + ^Fbihar (t) , 



(4) 



(5) 



with u = du/dr and F^ihar (t) = 77 sin (fir) + 
2 (1 — 77) sin (2r2r + (/j). Thus, Eq. (5) maps the pinion 
dynamics to the dynamics of a universal model - a Brow- 
nian particle moving on a periodic substrate subjected to 
a biharmonic excitation [23,2]. It is worth noting that, in 
spite of the abundance of numerical findings, the theoret- 
ical understanding of the directed transport phenomena 
represented by Eq. (5) remains far from being satisfac- 
tory [24] even about half a century after the earliest stud- 
ies [23]. The occurrence of DRT in Eq. (5) implies the 
breakage of two temporal symmetries: the shift symme- 
try and the time-reversal symmetry of the biharmonic 
excitation [2]. For deterministic ratchets subjected to 
biharmonic forces, it has been shown [18] that there ex- 
ists a universal force waveform which optimally enhances 
directed transport by symmetry breaking. Specifically, 
such a particular waveform has been shown to be unique 
for both temporal and spatial biharmonic forces. This 
universal waveform is a direct consequence of the degree- 
of-symmetry-breaking (DSB) mechanism. In particular, 
it is possible to consider a quantitative measure of the 
DSB on which the strength of directed transport by sym- 
metry breaking must depend. This mechanism has led to 
the unveiling of a criticality scenario for DRT. Indeed, it 
has been shown that optimal enhancement of DRT is 
achieved when maximally effective (i.e., critical) symme- 
try breaking occurs, which is in turn a consequence of 
two reshaping-induced competing effects - the increase 
of the DSB, and the decrease of the (normalized) max- 
imal transmitted impulse over a half-period - thus im- 
plying the existence of a particular force waveform which 
optimally enhances DRT (see [18] for additional details). 
Since thermal noise is significant in magnitude and un- 
avoidable at the nanoscale, the following question natu- 
rally arises: How does the DSB mechanism work at finite 
temperatures? 
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FIG. 2: (Color online) (a) Average velocity {{it)) = 
{{Adu/dr')) [cf. Eqs. (5) and (7)] versus relative amplitude 77 
for ip — ipopt = 7r/2, Q = O.OStt, 7 = 2, and three values of the 
noise intensity, (b) Normalized biharmonic function [Eq. (6)] 
versus time for Q = 1 and three values of rj. (c) Value of rj 
where the average velocity is maximum, rj'^pt^, versus a [cf. 
Eq. (5) with Fl,ii^^^ (r) instead of Fuhar (t), see the text] for 
1/5 = fopt = 7r/2, n = O.OStt, 7 = 2, and cr = 1. Also plotted 
is the theoretical prediction for the purely deterministic case 
Tjopt («) = 2a/ (1 -I- 2a) (dashed line) and the function 17* (a) 
(see the text, solid line). 



Here, we shall address this important question and 
provide analytical estimates for the dependence of the 
DRT on the system's parameters which are in excellent 
agreement with numerical results. To study numerically 
the effect of thermal noise (cr > 0) on the purely deter- 
ministic ratchet scenario, we calculated the mean veloc- 
ity on averaging over different realizations of noise {{u}} 
(cf. Eq. (5)). This corresponds to the mean velocity of 
the pinion ((i)) = A((m)) / (27r). Since Gaussian white 
noise does not break any relevant symmetry of Eq. (5), 
and the ratchet universality [18] predicts (for (7 = 0) 
that the optimal value of the relative amplitude rj comes 
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FIG. 3: (Color online) Average velocity {(u)) [of. Eq. (5)] 
versus relative amplitude r/ and noise intensity a for (p = 
ipopt = 7r/2, Q = O.OStt, and 7 = 2. The ranges of small 
and large noise intensities are shown in the top and bottom 
panels, respectively. 



from the condition that the amplitude of sin (fir) must 
be twice as large as that of sm{2ilT + ip) in Eq. (5), 
naively, one might expect that the average velocity of the 
pinion should present, as a function of 77, a single maxi- 
mum at ?7 = r/opt =4/5 when (p ^ 0, tt, as in the purely 
deterministic case. However, our numerical estimates of 
the 77 value at which the average velocity is maximum, 
77^^°, indicate a systematic deviation from ryopt = 4/5: 
A77 = 4/5 — r7o^° > 0, which is independent of the noise 
intensity over a significant finite range, as is shown in 
Figs. 2(a) and 3. To explain this paradox, the follow- 
ing remarks are in order. First, the effect of noise on 
the DRT strength depends on the amplitude of the bi- 
harmonic excitation while keeping the remaining param- 
eters constant. Second, while changing r\ and Lp allows 
one to control the breakage of the aforementioned rele- 
vant symmetries, it also changes the amplitude of the bi- 
harmonic excitation. Since the strength of any transport 
(whether or not induced by symmetry breaking) depends 
on the amplitude of the driving excitation, one concludes 
that these two effects of changing rj 01 ip overlap, so that 
one will find it difficult to distinguish the contribution to 
transport that is purely due to symmetry breaking, and 



hence to clarify the interplay between noise and sym- 
metry breaking. We shall therefore consider an affine 
transformation of the biharmonic excitation Fhihar (t), 
for the optimal value ip — Popt = 7r/2 for example [25], 
to change its image to [— 1/2, 1/2] , V77, thus making it 
possible to characterize the genuine effect of noise on the 
purely deterministic ratchet scenario: 
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8/9 < r/ ^ 1 (see Fig. 2(b)). After substituting Eq. (6) 
into Eq. (5), one straightforwardly obtains 
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= 7W^ + 7/:=V2(^') + V^C(t'), (7) 



where A = A (77) = M - m, r' = r' (r, 77) = At, = 
W [t]) = {M + rn) / {2k), Vt' = ^' {^^rf} = fl/A, and 
a' = a' {(J,ri) = a/A. Numerical experiments confirmed 
that the transport properties of Eq. (7) are similar to 
those of Eq. (5) in the sense that, for both equations, the 
average velocity {{u}) = {{Adu/dr')) presents a single 
extremum at the same value of 77 for a fixed set of the 
remaining parameters (see Fig. 2(a)). It is worth noting 
that the function A (77) is merely the width of the im- 
age of Fuhar [t), i.e., the difference between its maxima 
and minima as a function of 77, and that it presents a 
single minimum at 77 = 77* = 6/7. Also, the function 
W (77) represents an 77-dependent "load" force having a 
single maximum (in absolute value) at 77 = rjopt = 4/5, 
while (77 = 0, 1) = 0, as expected. Note that these 
particular values of 77* and 77opt are a direct consequence 
of the application of ratchet universality to the spe- 
cific form of the present biharmonic excitation F^ihar (t) ■ 
However, to better understand the roots of the present 
problem, it is convenient to consider the more general 
form F^.^^^ (r) = 77 sin (f^r) -|- a (1 - 77) sin (2^2t + ip), 
with a > being a parameter. For this case, one has 
77* = ri*{a),r]opt = 770^4 (a) ee 2a/ (1 + 2a), while the de- 
viation A7; (a) = rjopt (a) — r\%^t (o^) suggests a certain 
correlation between 77opt (a) and rfg^i [a) over a wide 
range of a values from a ~ 2 (see Fig. 2(c)). Together, 
these results therefore allow one to draw the following 
conclusions from Eq. (7). 

First, the aforementioned twofold transport effect of a 
biharmonic excitation F{,ihar (t), as 77 varies from to 1, 
may be decoupled into two terms: a constant excitation, 
W, and a biharmonic excitation, f*^^/2 (^')' having an 
amplitude which is independent of 77. The relevant obser- 
vation is that both excitations yield a maximum strength 
of transport at 77 = rjopt = 4/5. Therefore, replacing 
Fuhar (t) with /*^^/2 (t) m Eq. (5), 
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should yield a maximum average velocity at 7/ = rjopf = 
4/5, as is indeed confirmed by numerical experiments (see 
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Fig. 4). Thus, we propose for the system (8) the following 
scaling for the average velocity: 



(9) 



where C is a fitting constant that depends on the re- 
maining system parameters. Furthermore, for sufficiently 
high temperature (i.e., sufficiently far from the "steps" 
regime occurring at cr = 0, see Fig. 3) and driving ampli- 
tude (i.e., in the absence of stochastic resonance effects), 
exact agreement between numerical results and scaling 
(9) is expected over the complete range of r] values (see 
Fig. 4), while the scaling (9) remains valid over a wide 
range of frequencies (data not shown). 




FIG. 4: (Color online) Average velocity ((u)) [cf. Eq. (8)] 
versus relative amplitude rj for fl = O.OStt, "/ — 2, and two 
values of the noise intensity. Also plotted is the scaling law 
(9) for C = 1.00017 {a = 4) and C = 0.89072 (a = 0.5) [cf. 
Eq. (9); solid lines]. 

And second, alternatively to the case discussed in the 
first conclusion, only rescaling the temperature with the 
width A in Eq. (5) in accordance with Eq. (7), 



■it -I- sin M = vAct^ (t) + ^Fuhar (t) 



(10) 



should also yield a maximum average velocity sX rj = 
rjopt = 4/5. Numerical experiments also confirmed 
this prediction, as in the illustrative examples shown in 
Fig. 5(a). Thus, it is only after rescaling a — > Act in 
Eq. (5) that one recovers the purely deterministic ratchet 
scenario, which was an unanticipated result. To make the 
comparison between Eqs. (5) and (10) clearer, let us first 
transform Eq. (10) into the equation 



it + — sin u = a/ct£ (t) 
A 



+ 7' [77 sin (n'r) + 2 (1 - 77) sin (217't - 
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by rescaling the time, r — > At, and where 7' = 7/ A, Q' = 
0/A. Since the transport properties of Eq. (11) are sim- 
ilar to those of Eq. (10) in the sense that their respective 
average velocities present a single extremum at the same 
value of 77 for a fixed set of the remaining parameters, and 
that such an optimal value, rjopt = 4/5, is independent 
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FIG. 5: (Color online) (a) Average velocity (("")) [cf. 
Eq. (10)] versus relative amplitude rj for Q = O.OStt, 7 = 2, and 
three values of the noise intensity, (b) Value of rj where the 
average velocity is maximum, rj'^pt" (a), versus function rj*{a) 
[cf. Eq. (5) with Fl^^f^^^ (r) instead of Fbihar (t), see the text] 
and linear fit (12) (solid line) over the range 2.3 ^ a ^ 7.2 for 
V ~ fopt = 1^/2,0, = O.OStt, 7 = 2, and cr = 1. (c) Average 
velocity ((u)) [cf. Eq. (5)] versus initial phase difference ip 
for Q. = O.OStt, 7 = 2, cr = 1, and three values of rj. 



of the particular values of the amplitude and the fre- 
quency of the biharmonic excitation [18], one concludes 
from the comparison of Eqs. (5) and (11) that the ef- 
fect of thermal noise on the purely deterministic ratchet 
scenario can be understood as an effective noise- induced 
change of the potential barrier [d = d {rj) = A (77)) which 
is in turn controlled by the degree-of-symmetry-breaking 
mechanism through the function A (77) . Recall that the 
average velocity ((u)) exhibits, as a function of the po- 
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tcntial barrier d, a single maximum due to the thermal 
interwell activation (TIA) mechanism and the limiting 
behaviours lim£j_j.o,oo ((m)) = [24]. Obviously, the same 
scenario holds when Fbihar (t) is replaced with F^^-f^^^ (r) , 
which allows one to understand the behaviour of the de- 
viation At] (a) = Tjopt {oe) — riopt^ (o^) ^ ^ is changed (cf. 
Fig. 2(c)). Clearly, one can distinguish three regimes. 
Over the range < a < 1/2, the TIA mechanism dom- 
inates over the DSB mechanism. Indeed, the contribu- 
tion of the DSB mechanism to directed transport be- 
comes ever smaller as a ^ because ?7opt (a ^ 0) ^ 0, 
and hence the (whole) amplitude of the biharmonic ex- 
citation for which the DRT strength is maximum also 
becomes ever smaller as a 0, while the contribution 
of the TIA mechanism remains significant over the en- 
tire range < a < 1/2. Then one observes a transition 
regime over the range 1/2 < a < 2 as the effect of the 
DSB mechanism strengthens, which is manifest in the ex- 
istence of a narrow range of a values in which Ar] (q) « 
(see Fig. 2(c)). Finally, for the range a > 2, the contri- 
butions of the DSB and TIA mechanisms are correlated 
in the sense of the aforementioned effective noise-induced 
change of the potential barrier. This means that "rfg^i (a) 
is expected to be proportional to the value of the rela- 
tive amplitude where the effective noise-induced potential 
barrier presents a minimum: 

<>0(a)~^7?*(a)+B, (12) 

with A ~ 1.71, B ~ —0.744 being constants that are in- 



dependent of the remaining system parameters. The gen- 
eral scaling (12) is confirmed by numerical experiments 
(see Fig. 5(b)). Also, with regard to the dependence of 
the DRT strength on the initial phase difference, numer- 
ical results confirmed the scaling ((w)) ~ C" sin kp, where 
C" is a fitting constant that depends on the remaining 
system parameters, in accordance with ratchet univer- 
sality [18,19] (see Fig. 5(c)). 

In summary, we have explained the interplay between 

thermal noise and symmetry breaking in the ratchet 
transport of a Brownian particle moving on a periodic 
substrate subjected to a temporal biharmonic excita- 
tion, in coherence with the degree-of-symmetry-breaking 
mechanism. We have discussed how the present find- 
ings can be readily tested experimentally by means of a 
strongly damped rack-pinion ratchet at finite tempera- 
ture, where the rack undergoes a biharmonic lateral mo- 
tion. The general ratchet scenario presented in this Let- 
ter provides a theoretical framework for the optimal con- 
trol of the dynamics of nanomechanical ratchets in future 
applications. 
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